Two-photon lasing by a single quantum dot in a high-Q microcavity 
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We investigate theoretically two-photon processes in a microcavity containing one quantum dot 
in the strong coupling regime. The cavity mode can be tuned to resonantly drive the two-photon 
transition between the ground and the biexciton states, while the exciton states are far-off resonance 
due to the biexciton binding energy. We study the steady state of the quantum dot and cavity field 
in presence of a continuous incoherent pumping. We identify the regime where the system acts 
as two-photon emitter and discuss the feasibility and performance of realistic single quantum dot 
devices for two-photon lasing. 



I. INTRODUCTION 

Two-photon emission is an inherently nonclassical ef- 
fect. Realizing such emitters is important for instance 
in quantum telecommunication^ 2 -. In this respect, an in- 
coherently pumped source acting as two-photon laser is 
an attractive goal 3 -. Degenerate two-photon lasing, i.e., 
emission of two photons into the same mode of the elec- 
tromagnetic field, has been reported in atomic systems 
in two exemplary experiments, in the microwave range of 
the spectrum^ and later in the optical range, by pump- 
ing the dressed levels of potassium atoms 5 - - — . There is 
an extensive literature on degenerate two-photon lasing 
in atoms, some relevant references can be found in£~— . 
Further theoretical works discuss the perspectives of re- 
alizing two-photon lasing in semiconductor devices, such 
as in Refsii£~— . Interest from the solid-state community 
towards two-photon lasing has been renewed in view of 
recent experiments, demonstrating two-photon emission 
in semiconductor systems 2 ^— , and of the remarkable ex- 
perimental progress in achieving strong coupling between 
a single quantum dot and a microcavity mod o 23 ' 24 . In 
particular, it is expected that by tailoring the coupling 
between the quantum dot and the cavity, two-photon 
emission can be suitably enhanced, while one-photon pro- 
cesses are suppressed 3 -^. 

In this paper we theoretically analyze the possibility of 
observing two-photon lasing in a semiconductor device 
consisting of a single quantum dot strongly coupled to 
the mode of a microcavity. 

The discrete spectrum of a quantum dot (QD) supports 
a series of optical transitions that can be coupled to one 
or more cavity modes. Apart from a few localized state 
in the neutral QD, this nanostructurc can also accommo- 
date several carriers in several states, including the con- 
tinua of states in the wetting layer above the QD electron 
and hole barriers. When subject to Coulomb interac- 
tion, these electronic configurations will form a very rich 
N-exciton optical spectrum 2 ^—. In spite of this com- 
plicated structure, relevant for many properties, we will 



focus here on a "minimal" model in which just the low- 
est energy excitations of the neutral QD are considered 
as sketched in Fig. [1] where the QD is composed by four 
electronic levels: ground, excitons, and biexciton states 
(see Fig.QJb)). Previous theoretical works used this level 
scheme in order to study the emission of polarization- 
entangled photon pairs from a single quantum dot^Sr— . 
In the present paper, we consider a setup where the 
same cavity mode couples with the biexciton-exciton and 
exciton-ground state transitions, setting hence the basic 
requirements for the emission of two photons into the 
same mode of the electromagnetic field. 

Two-photon gain is reached in the regime where the 
resonator mode drives resonantly two-photon processes, 
while one-photon processes are far-off resonance. This 
is achieved by tuning the cavity frequency to half the 
biexciton-ground state transition frequency. Due to the 
binding energy of the biexciton state, this frequency 
can differ by several meV from the exciton transition, 
hence in principle offering the ideal conditions for re- 
alizing two-photon lasing^ 5 -. In contrast with previous 
studies, where population inversion is reached by ini- 
tially preparing the emitter in the excited state of the 
cascade emissio n 12 ' 13 ' 15 , here the quantum dot is con- 
tinuously pumped, and population inversion may result 
from the competition between coherent coupling, pump- 
ing and losses. 

Our study is based on a full quantum mechanical de- 
scription of the quantum dot and of the cavity field mode. 
We describe the dissipation caused by the cavity decay, 
radiative decay and pure dephasing of the quantum dot in 
the framework of a Born-Markov master equation. In ad- 
dition, we assume that the quantum dot is continuously 
pumped by an incoherent field. Incoherent continuous 
pumping is a typical excitation scheme of semiconduc- 
tors, where it has proved to be a determining ingredient 
in the description and achievement of strong coupling^ 5 -. 
The behaviour of the system is studied in the steady state 
of the dynamics, by analyzing the quantum-dot levels 
population, the cavity photon statistics, and the spectra 
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of the light intensity at the cavity output as a function 
of the various parameters. By means of this analysis 
wc identify optimal parameters for observing two-photon 
gain and we study the transition into the lasing regime. 
We carry out our analysis with realistic set of parameters, 
as compared to the semiconductor state of the art, and 
conclude that semiconductor devices may be suitable for 
two-photon lasing. 

This article is organized as follows. In Sec. HH we in- 
troduce the theoretical model. The coherent dynamics is 
summarized in Sec. Ill Al The incoherent processes are 
then introduced in Sec. Ill Bl The results are presented 
in Sec. IIII[ where we numerically study the steady state 
of the system, quantum dot and cavity field mode, as a 
function of the various parameters. Moreover, a semi- 
classical rate equation for the cavity mode, valid in the 
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FIG. 1: (Color online) (a) A high-Q microcavity coupled to 
a single quantum dot can be experimentally realized, for in- 
stance, by embedding the quantum dot in a photonic crys- 
tal. In (b) the ground (|G>), excitons \E 2 })), and 
biexciton (\B)) energy levels of the quantum dot are dis- 
played, the arrows indicate the coupling due to the resonator 
mode, the horizontal dashed line shows the frequency of the 
cavity-mode photon (see Sec.|lI]for details on the parameters). 
The cavity mode is resonantly coupled with the biexciton- 
ground state transition via two-photon processes. Due to the 
biexciton-binding energy shift, the exciton states are coupled 
far-off resonance. In this regime we explore the performance 
of the system as a two-photon laser. 



limit of a very good cavity and for the initial stages of the 
dynamics, is obtained in Sec. MI Dl The conclusions are 
drawn in Sec. IIVI where a detailed discussion on the ex- 
perimental feasibility of two-photon lasing with a single 
quantum dot in a cavity is presented. In the appendix 
the details of the calculations, at the basis of the rate 
equation in Sec. IIIIDI are reported. 



II. MODEL 

The system we consider is composed by a QD optically 
coupled to a single mode of the electromagnetic field of 
a microcavity with a high quality factor Q. The relevant 
states of the quantum dot are the four states of a biexci- 
ton system, sketched in Fig. [TJb): the ground state |G), 
at zero energy, the single exciton states \Ei) and I-E2) 
at frequencies oj\ and W2, respectively, and the biexciton 
state \B) at frequency 



(1) 



where x denotes the biexciton energy shift due to the 
exciton-exciton interaction. The transitions are coupled 
in a diamond-like configuration of levels with the sin- 
gle mode of the electromagnetic field, which couples to 
all dipolar transitions of the quantum dot. The master 
equation for the density matrix p of this system, QD and 
cavity-mode degrees of freedom, reads 



d 



(2) 



where H is the Hamiltonian describing the coherent evo- 
lution and L is the superopcrator which describes the 
effect of noise, dissipation, and incoherent pumping over 
the coupled system. The Hamiltonian H takes the form 

H = huj a a^a+ V fvjj j \Ej){Ej\+fvjj B \B){B\ 
3=1,2 

+hg J2 [(\G)(E j \ + \E j )(B\)a^+ll.c.] , (3) 
3=1,2 

with uj a the cavity frequency, a and the annihilation 
and creation operators of a cavity photon, and g is the 
vacuum Rabi frequency, which is here assumed to be 
equal for all transitions. It is important to note that we 
make this approximation in order to simplify the discus- 
sions in the following Sections. Typically, the couplings 
between each of the two excitons and the cavity mode 
differ, as they depend on the QD and cavity geometry 
(see, for instance, Ref. HH and references therein). How- 
ever, considering different coupling constants for the four 
transitions in the system is not significant for the two- 
photon lasing mechanism, as we show in Appendix B. 
There, we carried out a systematic analysis of this issue 
and we found that, within some reasonable restrictions 
for the parameters, the results are qualitatively the same 
as the ones we present in the main text using homoge- 
neous couplings for all the transitions. We also note that 
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this simplifying assumption is relevant for the coherent 
dynamics, as it introduces additional symmetries in the 
level couplings, see for instance Refs. I36ll37l 

As we are interested in two-photon processes, we will 
focus on configurations where the intermediate, exciton 
levels are far-off resonantly coupled to the radiation, 
while ground state and biexciton state are resonantly 
coupled by two-photon processes. In what follows we 
neglect the energy difference between the exciton states, 
assuming that the frequency difference \u>i — u 2 \ is much 
smaller than all the other characteristic frequencies of the 
system. In this limit, which is representative enough of 
more general situations, see App. B, we define the exci- 
ton frequency u>e = ~ ui2- In the frame rotating at 
the cavity frequency, the Hamiltonian becomes 



+hgJ2 [(|G)(^-| + |^)(B|)a+ + H. C .] ,(4) 
i=i,2 



where we have introduced the one-photon detuning 

2 ' 

and where 8 is the biexciton detuning, defined as 
8 = ujb — 2w Q ■ 



(5) 



(6) 



We remark that our model scheme assumes that the 
same cavity mode couples to both biexciton-exciton state 
transitions, \B) o- \Ej), and to both exciton-ground state 
transitions, \Ej) <H> \G). This choice is made in order 
to realize a degenerate two-photon laser, and requires a 
cavity where the two polarization modes have well sep- 
arated frequencies, as well as a proper definition of the 
quantization axis (which is naturally set by the geometric 
asymmetry of the quantum dot). 



A. Coherent two-photon processes 

We first focus on the coherent dynamics of the QD 
coupled with the microcavity field, discarding for the mo- 
ment any source of noise, and identify the region of pa- 
rameters in which two-photon absorption and emission 
are significantly enhanced over single-photon processes. 

We denote here by two-photon oscillations the coher- 
ent exchange of photons in pairs between the cavity mode 
and the quantum dot. This is found when excitations 
and de-excitations of the transition \G) O \B) corre- 
spond to resonant absorption and emission, respectively, 
of two photons of the cavity mode (corresponding to set- 
ting 2u) a — ujb, i.e. 8 ~ 0), while the single-exciton states 
remain effectively unoccupied. The latter condition is 
achieved when the single-photon detuning is much larger 
than the coupling to the cavity mode, x ^ 9 \J ( n ) > where 
(n) is the mean intracavity photon number. In this limit, 



we expect Rabi oscillations between the ground and the 
biexciton state with maximum amplitude, such that peri- 
odically one finds population inversion in the transition. 
The possible relevant cases are summarized in Fig. [5J 
where the cavity-mode frequency is shown with respect 
to the QD levels when the biexciton is bound. In par- 
ticular, the two-photon resonance regime is sketched in 
Fig. Oil) (see the caption for details). 

Energy/6 
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FIG. 2: (Color online) Energy levels of the QD for differ- 
ent cavity-mode frequencies, when the biexciton is bound 
(x > 0)- hi (I) the cavity mode is resonant with the tran- 
sition \G) <H> \E), namely, 8 — —\ (cue = u} a )- In (II) the 
\G) -o \B) transition couples resonantly with two photons of 
the cavity mode, 8 ~ (2ui a — IjJb)- In (III) the cavity mode 
is resonant with the exciton- biexciton transition \E) <H> \B) 
(8 = x, corresponding to uj a = 0Jb — oje)- 



Let us now describe in more detail the dynamics of 
this system at two-photon resonance. In order to pro- 
ceed with our analysis we introduce the basis of states 
|£, n), where £ = G, E\, E%, B denote the QD states and 
n is the number of cavity photons. Hamiltonian © cou- 
ples the closed set of states \G,n), \Ej,n—l), \B,n — 2). 
For 1 8 1 <C Q\J (n) <C x one can eliminate the exciton 
states from the coherent dynamics of the ground and 
biexciton states. By applying the Schrieffer-Wolf trans- 
formation^ one obtains an effective hamiltonian for the 
reduced Hilbert space of the QD states {|G),|B)} and 
for the cavity mode states, which reads 



h(S + 8 2P )\B)(B\ 

+hg cS fa 2 |i?)(G|+ a t 2 | G )( J B| 



(7) 



where 



'JcS 



x 



(8) 



scales the amplitude of the two-photon Rabi oscillations 
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between states \G) and \B), and 

4g 2 

S 2 p = — (9) 

X 

is the frequency shift on the two-photon transition, re- 
sulting from the coupling to the intermediate, exciton 
states. Correspondingly the ground-biexciton transition 
\G) <H> \B) and cavity field are in two- photon resonance 
when 

5=S 2P . (10) 

This frequency shift is independent of the number of pho- 
tons and has been discussed in the literature on two- 
photon lasers as self-induced frequency pulling terrn^. 
Hamiltonian (JT)) has been obtained assuming that x is 
the largest parameter, so that the intermediate states are 
coupled far-off resonance from ground and biexciton state 
by the cavity-mode photons^. We note that the ampli- 
tude of the two-photon coupling, Eq. ([8]) , is inversely pro- 
portional to x- Hence, by choosing large single-photon 
detunings, of the order of x, the effect of the interme- 
diate states is negligible. This choice implies, however, 
that the period of the two-photon Rabi oscillation, pro- 
portional to X: becomes very large, so that the effect 
of noise and dissipation, neglected in this analysis, may 
become relevant. 
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FIG. 3: (Color online) Populations of the QD levels p^{t) as 
a function of time, in units of 1/g, for the initial state \G, 2), 
obtained by evaluating numerically the equation p%(t) = 

I 1 2 

\(t;,n\e~ iHt/h \G,2)\ . The curves correspond tops (solid- 
blue line), pe = Pe ± + pe 2 (dotted-purple) and pg (dashed- 
brown). In (a) 8 — —20g, (b) 8 — ~2g, (c) 8 — —82P, while 
X = 20g in all three cases. See also Fig. [2] 



therefore the cavity frequency ui a ), corresponding to the 
situations sketched in Fig. [5] In Fig.[3Ja) the cavity is res- 
onant with the electronic transition between ground and 
single-exciton states, and we observe single-photon Rabi 
oscillations. In (b) the value of the cavity frequency cu a is 
far-off resonance from any electronic transition, and one 
observes that the population of the initial state is almost 
unit at all times. When the two-photon resonance con- 
dition (fit)]) is fulfilled, the coherent two-photon exchange 
\G,n) <-> \B,n — 2) is realized with maximum amplitude, 
as we can see in Fig. (He). Transfer of population be- 
tween ground and biexciton states occurs with a period 
of time of the order of t 2 p — n/(2g e g) (for n = 2), that 
corresponds to the time estimated from the dynamics 
predicted by Hamiltonian ([7]) , and found in perturbation 
theory by eliminating the intermediate levels. Note that 
the residual population of the single exciton state, which 
is observed in this plot, can be further suppressed by in- 
creasing the biexciton binding energy x, at the expenses 
of increasing t 2 p . The curve, corresponding to the ground 
state population, gives also the probability of finding two 
photons inside the resonator. 

Figure [4] displays the maximum value of the amplitude 
of the oscillation of the populations of the QD states 
as a function of S for x = 20_g, assuming that the initial 
state is \B, 0). Unit amplitude is found when the cavity is 
either resonant with the exciton-biexciton transition (5 = 
X, see Fig. [5] (III)) or with the ground-biexciton transition 
(S ~ 0, see Fig. [2] (II)). The widths of the two peaks are 
related to the strength of the excitons photon coupling. 
In particular, the two-photon peak is located at the two- 
photon resonance, Eq. (|10p . and is much narrower than 
the single-photon resonance because of the second-order 
coupling between ground and biexciton state. 



0.5 




-10 10 20 



6 (units of g) 



In Fig. |3] the time-dependent behavior of the popula- 
tion of the QD levels is evaluated numerically for differ- 
ent values of the cavity frequencies and assuming that 
at t = the system is in level \G, 2). The curves are 
found by integrating numerically Eq. @, setting the in- 
coherent term to zero. The biexciton binding energy is 
fixed at x = 20g. The Figure displays the time evolu- 
tion of the population of the ground, intermediate, and 
biexciton levels for different values of the detuning 5 (and 



FIG. 4: (Color online) Amplitude of oscillations of the popu- 
lation of the QD states, = max{p^(t)} — min{pf(t)} as a 
function of 8, in units of g, with \ — 20g and initial state 
\B,0). The curves correspond to Vb (solid-blue line), Ve 
(dotted-purple) and Vg (dashed-brown). 
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B. Incoherent processes 

We now include the incoherent processes in the dy- 
namics of QD and cavity mode. In the present model, 
they are due to the coupling of system with the environ- 
ment constituted by the external modes of the electro- 
magnetic field and the phononic modes of the sample in 
which the QD is embedded, and are described by the Li- 
ouvillian L in the master equation, Eq. ([3]). This includes 
the decay of the cavity-mode photons into the external 
modes of the electromagnetic field at rate k, decay of the 
exciton and biexciton states at rate 7 and 75, respec- 
tively (which include radiative and non-radiative losses, 
the latter due to electron-phonon interactions), and pure 
exciton dephasing at rate T^, due to interaction with the 
bulk phonons, which we assume here to be MarkovianiS. 
In addition, the QD is pumped continuously. Typically, 
electron-hole pairs are created in high energy states by 
electrical injection or optical pumping. Then they relax 
to the excitonic levels of the QD producing an effective 
incoherent source of particles with injection rate which 
we denote here by P. The Lindblad term L in Eq. (|3|), 
describing these processes, takes the form 



■r^-C(\E s )(B\) 



3 

D, 



53^(l^i>(^il)p + r*c(|B>(B|)p 



--^[CQEjXGD + CdBXEj^p, (11) 



where we have used the definition 



C{x)p = 2xpx' 



x'xp 



px' x , 



(12) 



with x an operator acting on the Hilbert space of cavity 
mode and QD. In the following we assume for simplicity 
that 



IB = 27 . 



(13) 



We finally note that in the literature pure dephasing 
has been sometimes described with non-Markovian type 
of Liouvillian o 37 ' 41 ' 42 . It is indeed expected that non- 
Markovian corrections will give rise to non-linearities in 
the phonons-cavity coupling. In this paper we describe 
dephasing as Markovian process 4 ^, as our scope is to 
identify the parameter regime where two-photon lasing 
can be observed. 

The incoherent dynamics in general competes with 
the coupling with the cavity mode. Noise and dissipa- 
tion, hence, will tend to suppress two-photon coherent 
exchange between QD and cavity mode. Nevertheless, 
quantum effects, at the levels of few photons, will af- 
fect relevantly the dynamics provided that the coupling 
strength of the QD with the cavity mode exceeds the loss 



rates. An important parameter, determining whether the 
system is in the strong coupling regime, is the coopera- 
tivity. It is defined a o 43 ' 44 : 



C 



2 jL 

Ky 



(14) 



where 7' = 7 + IV In absence of incoherent pumping 
strong coupling is warranted when C > 1. 



III. TWO-PHOTON LASING 

In this section we study the steady-state properties of 
the system evaluated by solving numerically the master 
Eq. ((2J) with the Lindblad term (fTTj) as a function of the 
physical parameters. We first focus on the case in which 
dephasing is set to zero, in order to single out the role 
of the various incoherent processes on the lasing dynam- 
ics, and then consider the effect of dephasing within the 
description in Eq. (fTTj) . We finally discuss a specific pa- 
rameter regime where an analytic solution can be found 
for the dynamics in a very good cavity, and allows us 
to identify a condition when two-photon gain exceeds 
one-photon gain in the initial stages of the dynamics. 
Comparison with the numerical results shows that this 
condition does not necessarily imply that the system ac- 
cesses the two-photon lasing regime. 

We remark that previous theoretical works studying 
two-photon lasing in cavities focussed on setups where 
atoms in the upper excited state were pumped into the 
cavity. In that case, population inversion on the two- 
photon transition is warranted by the state preparation of 
the atoms entering the cavity, and the rate at which they 
cross the resonator is one important parameter control- 
ling whether two-photon lasing is foun d 12 i 13 ' 15 . In this 
work we consider a single dot inside a resonator whose 
levels are continuously pumped by external, incoherent 
fields. Population inversion, if found, hence results from 
the competition between the pump, the loss and the emis- 
sion processes. 



A. Steady state and photon statistics 

We first study the steady-state population of the QD- 
electronic states as a function of the cavity frequency and 
of the pump strength. As is visible in Fig. O the cavity- 
QD system possesses three resonances: the two-photon 
resonance at 8 = —82P, when the transition \B) -H> \G) 
is resonantly driven, and the one-photon resonances at 
S = x an d 8 = — x, when the transitions \B) \Ej) and 
\Ej) -f-> |G), respectively, are resonantly driven. 

When the cavity is far-off resonance from these tran- 
sitions, then the steady state of the electronic levels 
is determined by the incoherent processes of exciton 
pumping and decay described in Eq. (|11[) . Denoting 
by pj — (j\p\j) the population of the electronic state 
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FIG. 5: (Color online) Steady state populations of the QD 
energy levels as a function of the cavity detuning S (in units 
of g) for n = 0.1<?, Td = 0, 7 = 0.05g, \ — 20g and pump 
rate (a) P = 0.1%, (b) P = 0.5ff, (c) P = O.85, (d) P = I.85. 
The curves correspond to the occupation of the states \B) 
(solid line), \Ej) (dotted line), and \G) (dashed line). Two- 
photon resonance is at S = — 82P ~ 0. At 8 = 20g transition 
\B) I Ej) couples resonantly with the cavity mode. 



FIG. 6: (Color online) (a) Average photon number (n), (b) 
Fano factor F, and (c) g (2) (0) as a function of the cavity 
detuning 5 (in units of g) for « = O.lg, Td = 0, 7 = 0.05g, 
X = 20g and pump strengths P — 0.16g (dashed lines), P = 
0.5g (solid lines), P = 0.8g (dotted lines), P = 1.8g (dotted- 
dashed lines). The horizontal dotted lines indicate (b) the 
value F = 1, (c) the values g^fO) = 1 (Poisson distribution) 
and #' 2 '(0) = 2 (thermal distribution). 



\j) (j = G, B, Ei, E2), in this limit they take the form 
p G = 7 2 /(F + 7 ) 2 , PB = P 2 /{P + 7 ) 2 , and p El = Pe 2 = 
jP/(P + 7) 2 , showing that the biexciton state is occu- 
pied with almost unit probability when P exceeds the 
radiative decay rate. At the resonances, this behaviour 
is significantly modified by the competition between the 
incoherent processes and the strong coherent coupling to 
the cavity mode. 

Figure [5] displays the value of the populations at steady 
state as a function of the two-photon detuning S, where 
here d = ±20<7 corresponds to the one-photon resonances. 
The subplots from top to bottom are calculated for in- 
creasing values of the pumping rate P. Here, one ob- 
serves depletion of the biexciton population at S ~ 0,x, 
which is more enhanced at two-photon resonance. At 
8 ~ when the cavity mode drives resonantly the 
biexciton-exciton transition, depletion of the biexciton 
state population is accompanied by a significant increase 
of the intermediate levels occupations, while the ground 
state remains essentially empty. Photon emission on the 



biexciton-exciton transition is here expected. At two- 
photon resonance one observes a significant population 
of the ground state, while the (small) increase of the in- 
termediate states populations further indicates that the 
dominating processes are emission of photons in pairs 
from the biexciton to the ground state level. Increase 
of the pump strength leads to stationary population in- 
version between the biexciton and the other levels, and 
to a broadening of the interval of values of 6 in which 
population depletion is observed. 

We now analyze the cavity field properties as a function 
of 5. Figure [6] displays (a) the mean number of photons, 
(b) the Fano factor, which is the ratio of the variance of 
what is measured to that for a Poisson distribution with 
the same mean photon number, and it is defined as^ 



with n = a T a the photon number operator, and (c) the 
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FIG. 7: (Color online) Same as in Fig. [6j but for a fixed value 
of the pump rate, P — 0.5g, and for different values of the 
radiative decay rate: 7 = 0.05g (solid line), 7 = O.lg (dot- 
dashed line), 7 = 0.5;; (dashed line). 



FIG. 8: (Color online) (a) Average photon number, (b) Fano 
factor, and (c) g 2 (0) as a function of the pumping rate P 
(in units of g) for 5 = 0, Td — 0, 7 = 0.05.g, \ = 2Qg, and 
different values of the cavity decay rate: k = 0.2g (dashed 
line), k = O.lg (solid line), and k = 0.05g (dot-dashed line). 



intensity-intensity correlation function at zero delay, 

^(0) = ^=^, (16) 
(n) 

quantifying the probability of two-photon coincidence de- 
tection. The curves are calculated at relatively small 
values of the radiative decay and zero dephasing rate, 
and for different values of the pump strength. The pho- 
ton number exhibits two peaks, at two-photon and one- 
photon resonance, whereby the maximum at the two- 
photon resonance is higher than at one-photon resonance. 
This behavior is enhanced as the pump strength is in- 
creased. Moreover, while in general the peak at two- 
photon resonance is narrower than the peak at one- 
photon resonance, it broadens significantly for large 
pump powers. The photon-number fluctuations, on the 
other hand, decrease at two- and one-photon resonance. 
Fano factor and <7^(0) approach unity at 8 = x f° r 
large pump rates, while they remain larger than unity at 
5 = 0. Close to the minimum at two-photon resonance 
large fluctuations in the photon number are observed. Fi- 
nally, far from the one- and two-photon resonances the 
second-order correlation functions approaches the value 



g( 2 )(0) ~ 2, which indicates a thermal character of the 
field. 

The effect of radiative decay on the properties of 
the cavity field is analyzed in Fig. [7J Here, we ob- 
serve that when the decay rate is increased, the number 
of photons at one- and two-photon resonance decreases 
(while the photon number at the one-photon resonance of 
the exciton-ground state transition increases), and cor- 
respondingly the photon-number fluctuations increase. 
When P ~ 7, in particular, the Fano factor is maxi- 
mum at the two- photon resonance, see Fig. EJb). Corre- 
spondingly, g( 2 '(Q) ~ 2, indicating that the cavity field 
approaches a thermal distribution. 

A further important parameter is the rate at which the 
cavity field leaks out of the resonator. In Fig. [8] we study 
the behavior of the steady state photon statistics as a 
function of the exciton pumping rate for various values 
of the field decay rate. While it is obvious that, increas- 
ing the pump and decreasing the cavity loss rate, the 
number of photons inside the cavity increases, the be- 
haviour of the photon-number fluctuations is interesting. 
We observe that </ ' > 1 (F > 1) for all values of k 
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photon statistics, showing that the one-photon resonance 
as a function of 5 is higher and broader for the largest 
value of the dephasing rate here considered, Td = 0.5g. 
We note, however, that this effect is observed only when 
the dephasing rate exceeds the radiative decay rate, as 
we assume in this specific case. 
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FIG. 9: (Color online) (a) Cavity-field average photon num- 
ber (h) at steady state as a function of the dephasing rate Td 
in units of g. The solid line corresponds to S — (two-photon 
resonance), and the dashed line to S = \ (one-photon reso- 
nance), (b) Average photon number, (c) Fano factor, and (d) 
g^ 2 \0) as a function of the two-photon detuning S in units 
of g. The dot-dashed line corresponds to Td = 0, the solid 
line to Td — O.lg, the dashed line to Td = 0.5g. The other 
parameters are P — 0.5g, 7 = 0.05g and k = O.lg. 



and of P here considered, indicating that the field is su- 
pcrpoissonian. Moreover, the Fano factor increases with 
increasing pump rates at n = 0.2g, while it decreases to 
values which approach unity for k = 0.05(7 an d O.lg. ft 
is interesting to note that for a certain interval of values 
of P, which are much smaller than the cavity coupling 
strength g, the Fano factor evaluated at k = 0.05g and 
k = O.lg exhibits a broad maximum. We note that val- 
ues of the Fano factor F > 1 are expected in two-photon 
lasers, see for instance Re fi 13 ' 15 . A direct comparison 
with the values found in the literature and the ones we 
obtain is not trivial, as differing from the systems i n 13 ' 15 
the setup we consider has relatively large loss rates k and 
a different pumping mechanism. 

The effect of dephasing on the state of the cavity field 
is studied in Fig. [HI Figure [SJa) displays the mean intra- 
cavity photon number as a function of the dephasing rate, 
showing that (n) decreases at two-photon resonance as 
is increased, while it increases with Yd at one-photon res- 
onance. Figures Elb)-(d) analyzes this behaviour in the 



B. Spectrum of the cavity field 

We now focus on the spectrum of the field at the cavity 
output, which is evaluated from the formula 



7r(ata) c 



-ft 



(at(r)a(0)> c 



r dr . 



(17) 



using the quantum regression theorem^—, and where 
we assumed that the input field is in the vacuum, con- 
sistently with master equation @. We define = 
Trj-poo}, where p x is the density matrix of QD and cav- 
ity field at steady state. Below we analyze the spectrum 
when the system is driven at 8 — 0. 

Figure [TU] displays the spectrum at large loss rates and 
small cooperativities. We note that it exhibits three main 
peaks at the bare system frequencies: one at the cavity 
frequency w — w a , a second at the ground-exciton tran- 
sition frequency and a third one at the frequency of the 
exciton-biexciton transition. In particular, the central 
curves are Lorentzian with linewidth for the cavity 
field in absence of the QD, see Fig. [10] (c). In this regime, 
hence, the spectrum provides a spectroscopic measure of 
the characteristic frequencies of the system and of their 
corresponding lincwidths. We note, however, that due to 
the small coupling the exciton resonances are two orders 
of magnitude smaller than the central peak. The effect 
of increasing the pump power, nevertheless, is to increase 
the number of cavity photons at the frequency of the 
biexciton-exciton transition, while the number of pho- 
tons at the frequency of the exciton-ground state transi- 
tion becomes smaller. 

Significant differences from this behavior are observed 
at large cooperativity, with coupling strengths overcom- 
ing the loss rates. Figure [TTJ displays the spectrum evalu- 
ated as a function of the pump power and large coopera- 
tivities. The central peak is dominant and, by increasing 
the pump, moves from being centered at uj = uj a — §2P 
towards u> = ui a . As the pump power is increased, more- 
over, the exciton resonances both decrease, while the 
linewidth of emission line at the central peak becomes 
narrower, as it is visible in Fig. [TT] (c) , thereby exhibiting 
a typical signature of lasing at this frequency. 



C. Non-classical properties of the cavity field 

The analysis done so far shows several features that can 
be interpreted altogether so that the system accesses two- 
photon lasing at large cooperativity and sufficiently large 
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FIG. 10: (a) Spectrum of the cavity field intensity at the output for different values of the pump rate P with n = lOg, 5 = 0, 
X = 2000(?, 7 = 5g. The curves refer to P = g (dotted line), P = 5g (dotted-dashed line), P — 50g (dashed line), P — lOOg 
(solid line). The plot in (b) shows the curves in the interval of frequency close to tJ — tJ a = and plotted in linear scales. 
The thick gray curves in (a) and (b) indicate the (Lorentzian) emission line of the cavity in absence of the coupling with the 
quantum dot. The plot (c) shows the linewidth at half height of the peaks in (b) as a function of P. The horizontal dashed line 
indicates the value of k. 
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FIG. 11: Same as Fig.QJ)] but with k = 0.1#, S = 0, x = 20g, 7 = 0.05ff and P = O.Olp (dotted line), P = 0.05p (dotted-dashed 
line), P = 0.5<? (dashed line), P = 1.8g (solid line). 



pump power. In fact, at two-photon resonance and suffi- 
ciently large pumping rate we observe population inver- 
sion on the biexciton-ground state transition (as well as 
on the biexciton-exciton transition). Moreover, at these 
values the statistics of the light is superpoissonian, ex- 
hibiting values of the Fano factor which are consistent 
with the ones expected for a two-photon lase r 13 i 15 . In 
addition, for these parameters we observe that the cavity 
spectrum has intense emission and line narrowing at the 
cavity frequency. 

Two-photon lasers are expected to be nonclassical light 
sources, however, the analysis presented so far does not 
show clear nonclassical features of the cavity field. In 
the literature, nonclassical features of degenerate two- 



photon lasers are usually associated with squeezing of a 
quadrature at the cavity output. When the system is 
pumped by atoms prepared in the excited state, squeez- 
ing of the cavity field can be observed if a seed is put 
in the system. If, for instance, the cavity field has been 
initially prepared in a coherent state, squeezed light is 
observed when incoherently pumping the system for a 
transient time, of the order of typical time scale deter- 
mining phase diffusions^. In this work we check that 
such behavior is observed also in the setup we consider, 
namely, a single QD continuously pumped and strongly 
coupled to the cavity mode. The results are obtained by 
numerical integration of the initial master equation (jlip . 

Figure [redisplays the time evolution of the variance of 
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FIG. 12: (a) Time evolution of the variance (AX(<f>) = 
(X 2 ) — (X) 2 ) of two orthogonal quadrature of the field, de- 
fined as X = ae^'^+a' e 1 ^ (the horizontal dotted line indicates 
the shot-noise limit). The solid (dashed) lines corresponds to 
the quadrature AX((/>)| m i n (AX(0 + 7r/2)| max ) with minimum 
(maximum) value at the given time. The corresponding phase 
4> of the quadrature with minimum variance is plotted in (b). 
The field is initially prepared in a coherent state with aver- 
age photon number (n) = 4. The black lines are found for 
P = O.Olp, 7 = O.Olg and k = 0.01, whereas the gray lines 
are found for larger rates of incoherent processes: P = O.lg, 
7 = 0.05g and k = 0.02. The other parameters are x = 50g, 

S = — Szp. 



the field quadrature which possesses minimum variance 
at all times. For sufficiently short times, the field be- 
comes squeezed. At longer time the squeezing disappears 
because of the phase diffusion induced by the incoherent 
processes. The interval of time in which the squeezing is 
visible is shorter for larger rates of the incoherent pro- 
cesses. 



D. A rate equation for the cavity field dynamics 

The possibility of realizing media exhibiting two- 
photon gain have been extensively discussed in the lit- 
erature. Theoretically, one typically estimates gain as- 
sociated with two-photon processes using a semiclassical 
analysis. Two-photon lasing is expected when the two- 
photon gain exceeds the one-photon gaini^. 

In this section we derive a rate equation for the cavity- 
photon state occupations. The rate equation is obtained 



assuming that the intermediate (exciton) levels are cou- 
pled far-off resonance from the ground and the biexciton 
levels, it is valid for the initial evolution, when the intra- 
cavity photon number does not saturate the QD transi- 
tions, and takes the form 



Pn = 



- 1 ™ Pn 



G 



(2) 



2 Pn- 2 



+1 n+ lPn+l + 1 n+2 Pn+2 , 



(18) 



where p n denotes the probability that the cavity mode 
contains n photons. The details of the derivation are 
given in Appendix |XJ The rates take the form 

G« = (n + l)Ai , 

= (n + l)(n + 2)A a , 

r« = nC 1 , 

r« = n(n - 1)C 2 , 

r-^GW + G^+rW+if) 



(19) 
(20) 
(21) 
(22) 
(23) 



and they describe transitions connecting different pho- 
tonic states, as shown in Fig. [13] The coefficients ap- 
pearing in the expressions are rates. In particular, A\ 
and A 2 are the one- and two-photon gain, respectively, 



Ai = 



Ao 



8g 2 P [P(P + 3T d ) + 7 (6P + r d ) + 7 2 



Z2g A P 2 



\ 2 2 

•7) X 



(P 



, (24) 

iY{p + i + v d ) x 2 

while C\ and C 2 are rates of dissipation of one and two 
photons, respectively, and read 



G 2 



P' 



Ao 



(25) 



Equation (|T8f has the form of the standard rate equa- 
tion for lasing^, where in addition it includes two-photon 
processes. Differing from previous treatments on two- 
photon lasing, based on injecting atoms into a resonator 
prepared in the upper states, rate equation (|18p is de- 
rived for a medium which is continuously pumped inside 
the resonator. This approximation is adequate when the 
biexciton transition is not saturated, i.e. for low exci- 
tation of the system. In particular, the coefficients of 
the one-photon processes correspond to the coefficients 
of the laser rate equation (see, for instance, Ref. 1491 ) here 
expanded to lowest order in the field-medium coupling. 

Comparison between the different rates provides in- 
sight in the initial dynamics of the cavity field. In partic- 
ular, gain overcomes dissipation when the pumping rate 
exceeds the losses, P > 7 , n (we recall that this rate 
equation has been obtained in the limit of small K, see 
Appendix Moreover, two-photon gain will be (ini- 
tially) dominant over one-photon gain when the ratio 



G 



(i) 



(n + 2)A 2 
A x 



> 1 



(26) 
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FIG. 13: Sketch of the processes considered in the rate equa- 
tion (|18[) . The horizontal lines indicate the energy levels of 
the cavity mode, the vertical arrows indicate the flow of prob- 
ability in and out of the energy level at n photons with the 
corresponding rate. 



For n = 0, this reduces to the inequality 

2A 2 8g 2 



1 < 



-4, 



p2 



(27) 



which has been obtained in the limit P 3> 7,1^, show- 
ing that two-photon gain overcomes one-photon gain pro- 
vided that the cavity-exciton coupling strength is larger 
than the pumping rate. We remark that condition (|27j) 
has been obtained assuming a very good resonator, in 
the limit in which the cavity loss rate is the smallest 
parameter of the dynamics. We further note that con- 
dition (|2"T)) is consistent with the semi-classical result in 
Ref obtained for a quantum well, when studying two- 
photon amplification of a weak classical probe. Com- 
parison with the numerical results at large cooperativity, 
however, shows that the signatures of two-photon lasing 
are found when the pumping rate exceeds the coupling 
strength, contrarily to condition (|27|) . We attribute this 
discrepancy to the fact that lasing in a high-Q cavity is 
thrcsholdlcss 5 -^, therefore, the comparison of the two- and 
one-photon gain at the initial stage of the dynamics docs 
not provide sufficient information about the behavior at 
steady state. 



IV. DISCUSSION AND CONCLUSIONS 

We have analyzed the conditions for two-photon gain 
and lasing with a single quantum dot inside a high-Q 
microcavity. We have shown that, by taking advantage 
of the biexciton binding energy, the system can exhibit 
clear features of a two-photon emitter when the cavity 
frequency is tuned to the two-photon resonance. Two- 
photon lasing, in particular, is found in high-Q cavity for 
sufficiently large pump rates. 



A crucial point to be stressed is the feasibility of two- 
photon lasing by using QD embedded in semiconductor 
microcavities. Different techniques have been developed 
for fabricating systems in which just a single QD is lo- 
cated close to the maximum (antinode) of the electro- 
magnetic mode of a cavity built up either in a photonic 
crysta l 23 i 38 ' 51 ' 52 or in a micropillar— ~— . In order to es- 
tablish whether or not we have been working within a 
reasonable regime of parameters, wc will compare our 
values with the ones used in two of the most recent 
experiment: ] 38 ' 52 with incoherent pumping as it is the 
case of our theory. One of the strongest coupling re- 
ported up to now^ is g = 150/icV in a system with 
k = 83/icV = 0.55(7 an d 8 in the range from 100 to 
200/icV i.e. 5 « g. Smaller but similar values are those of 
Ref^ where g = 68^eV, k = 39^eV = 0.57g. These sets 
of parameters are not far from the ones we use for most 
of the results we present here which have been obtained 
for a ratio n/g = 0.1 which corresponds to 17% — 18% 
of that achieved in the experiments. Devices with such 
characteristics are expected to be realizable in the near 
future 5 ^. 

Taking 7 < k (see figures [5] or ITT|) is also very reason- 
able compared with currently existing systems. A cru- 
cial magnitude for getting two photon lasing seems to 
be the biexciton binding energy x that we typically take 
as 20 times g. With the experimental values above men- 
tioned for g this means x between 1.4meV and 3 meV, 
which are very typical values for self-organized growth 

qds. 

The comparison of the pumping rate P with exper- 
imental values is more complicated. One must con- 
sider two issues: the evaluation of the externally ap- 
plied power excitation and its reduction to the absorbed 
pumping. The first factor, the externally applied pump 
Pext, is just the number of photons arising to the sample 
per unit time. This can be evaluated as the pumped 
energy per unit time, i.e. the pump power, over the 
energy of one of those photons. In the case of Rcf^S 
a pump power of 1 nW and a wavelength of 818nm 
gives P ex t = 0.041ps _1 which in units of energy becomes 
P ext = l70fj.eV. In§S the wavelength is 800 nm while 
pump power is varied from 25 nW to 500 nW, which 
gives 418/icV < P ex t < 10474/icV. In the second step, 
getting P from P ext requires to estimate the reduction 
due to the reflection of the light by the sample, the ab- 
sorption losses in mirrors or the part of the light which 
is actually absorbed by the cavity itself. This reduction 
from external to absorbed pumping must be evaluated 
for each particular case but a reasonable estimation for 
the reduction is to get P between 1% and 5% of Pea t 58 ' 59 - 
This means that in the case of reference 3 ^ one can esti- 
mate P w 4.5/ieV = 0.03<7 while in the case of reference 52 
12/j.eV < P < 300/icV which means O.I8.9 < P < 4.5g. 
The results we present in this paper are obtained with 
0.0 lg < P < 5g which are values precisely in the regime 
of the experiments discussed above. 

The conclusion from the previous analysis shows that 
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the two photon lasing, as depicted by our results, can be 
feasible with currently available samples and technology 
for QD embedded in microcavitics. 
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Appendix A: Effective Dynamics and rate equation 

In this appendix, we first derive the effective master 
equation for QD and cavity mode in the regime in which 
X is the largest parameter. We start from the master 
equation (JTTJ) , and use the basis of symmetric and anti- 
symmetric superposition of the exciton states, 



|±) = i=(|£ 1 )±|£ 2 )), 



(Al) 



noting that the cavity mode couples with the symmetric 
state. In this basis Hamiltonian ((4]) takes the form 



obtaining the set of equations 

Pgg = -2PPGG + J(P+++ P — ) 

+ - 03 + iffeff) [a^apcG + a t_ Pbg - atp ++ a\ 
+H.c] + ^C(a)p GG , 

pBB = -2jp B B + P(P++ + P ) 

+ [- (a + H7cff) (aa^pBB + a 2 p G B - ap ++ <J) 
+H.c] + ^C(a)p BB , 

PGB = -(P + l + T d -i5)pGB 

- (a + ig e g) (atapGB + pbb - aV++ flt ) 
-ifi- iffcff) (pgb<i<J + PGGa f 



a) p ++ a) 



Pbg 
P— 

P++ = 



+ -C{o)pgb , 
Pgb > 



P+-) 



+PPGG + IPBB + Y P++ + ^ C ^P-- ' 

P + 7 + y) P++ - i|(P-+ - P+-) 

+PPGG + IPBB + -yP 

+ [—(/?- iffcff) (aoV++ _ apGBa - apGGa*) 
- {a - iffeff) (a)ap ++ - a) p B G^ - a^p B Ba) 
+H.c] + %C{a)p++ , 

+2 (/3 + iff eff ) p_ + aa t + 2 (a + ig eS ) p- + a)a 



-i|(P— - P++) + 2 £ ( a )^-+ 



P+- = P-+ 



(A3) 



where pij is here an operator over the Hilbert space of 
the cavity mode, and 



H 



X + S 



(I+X- 



-)(-{) + S\B)(B\ 



+v / 2 5 [(|G)<+| + |+)(i?|)at+H.c.] , (A2) 



where e = uj\ — u)2 is the frequency difference between 
the exciton states, which we assume |e| <C X- I n the limit 
in which x 1S the largest parameter, the single-photon 
coherences pj.G(n,m) = (j, n\p\G, m), and pB,j(n,m) = 
(B,n\p\j,m) (with m — n, n + 1) oscillate very rapidly 
with respect to the other elements of the density matrix, 
and can be adiabatically eliminated from the dynamics, 



3cff = -4— 



a 2 (P + 3 7 + 3r d ) 



g 2 (3P + 7 + r d ) 

9 

X 



(A4) 
(A5) 
(A6) 



The dynamics these equations describe is valid up 
to second order in the expansion in l/|xl- Moreover, 
we assume that the cavity decay rate k and the fre- 
quency difference |e| are the smallest parameters, such 
that k, |e| < g e s . Hence, the regime of validity of these 
equations is quite small and experimentally demanding. 
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FIG. 14: (Color online) Average cavity photon number (h), 
Fano factor F and <? 2 (0), evaluated numerically from Master 
equation (|ll|l (solid lines) and analytically using rate Eq. (|18[) 
(dashed lines) as a function of the pumping rate P (in units of 
g). The parameters are \ = 50t;, S = <$2P, = 0, 7 = 4p and 
k = 0.02<;. When the pump is sufficiently small, P < 7, such 
that the number of intracavity photon does not saturate the 
biexciton-ground state transition, the rate equation provides 
a good description of the dynamics at all times. In this regime 
the system behaves as a laser below threshold, and the cavity 
field approaches a thermal state. 



Nevertheless, they lead us to an analytic solution, which 
give us some insight into the process of two-photon las- 
ing. 

In the following we assume e = 0, so that the last two 
equations in Eqs. (|A3|) are decoupled from the rest, and 
the system dynamics can be described by the first six 
equations of Eqs. (|A3|) . We first note that the coher- 
ent dynamics described by this set of equation, obtained 
by setting P = 7 = r f / = n = 0, is equivalent to the 
two-photon dynamics described by the effective Hamilto- 
nian derived in Eq. ([7]). When considering the incoherent 
processes, on the other hand, one finds that the exciton 
states are involved in the system dynamics in a nontrivial 
way. 

We now aim at getting further insight into the regime, 
in which gain can be observed, and in particular in which 
two-photon gain is dominant. At this purpose we derive 
from Eqs. (|A3j) a rate equation for the photonic state 
occupations. We assume P + 7 ^> \g c s\. In this limit the 
atomic population and the two-photon coherences decay 
very fast with respect to the typical time scale of the 
evolution of the cavity field. Denoting by 

p n = Tr{\n){n\p} (A7) 
the population of the photon state |n), one obtains 



the rate equations for the population p„, reported in 
Eq. (fTgj) by expanding Eqs. (|A3|) at lowest relevant order 
in g c g/(P + 7) and taking the corresponding trace. 

We remark that rate equation (fTB)) has been obtained 
assuming very small cavity decay (high-Q resonators) by 
an expansion in the QD-cavity field interaction strength, 
and therefore it does not contain the higher-order, non- 
linear terms which describe the lasing regime. It is valid 
when these terms are negligible, namely, either (i) at all 
times when the pump is small such that the system be- 
haves as a laser below threshold, in which the field is 
thermal (see Fig. HI)) , or (ii) at the initial times, when 
the number of cavity photons is still sufficiently small 
that the nonlinear terms are negligible. 



Appendix B: Validity of the model 

In this appendix we study the limit of validity of the 
model introduced in section [TXJ In particular we focus 
on the simplifying assumptions of equal coupling con- 
stants between the cavity field and the four electronic 
transitions and on the assumption of negligible detuning 
between the single exciton states. 

In general the exciton states energies are not degener- 
ate. We denote their detuning by e. Moreover, we denote 
the coupling strengths between the cavity mode and the 
four electronic transitions by the symbols gci, Qgi-, 9ib 
and giB- The Hamiltonian, describing the coupling be- 
tween QD and cavity field, takes the form 

H = ^[(x + S + e)\E 1 )(E 1 \ + ( X + 5-e)\E 2 )(E 2 \} 
+h5\B)(B\ 

+hg G1 [(\G) (E 1 \ + a\G) (E 2 \) J + H.c] 

+hg 1B [(\E 1 )(B\ + p\E 2 )(B\)ai +H.c] , (Bl) 

where we have introduced the dimcnsionlcss coefficients 
a and /3, giving the ratios a — gcij gc\ and /3 = g 2 s/ 'giB- 
Under the assumption of large x^ we eliminate the inter- 
mediate, exciton states and obtain the effective Hamilto- 



H eS = h[6 + S2p + 5 n a J! aj \B)(B\ 
+g cfi (a 2 \B)(G\+a^\G)(B\ 



(B2) 



with 



5 2 p 
5cff 



(1 + /3 2 ) 



X 



ZgaigiB 

X 



9ib0- 



1 



(B3) 



Differing from the ideal case discussed in this article, 
where we take all coupling strength to be equal, we ob- 
serve that in general the light shift of the two photon 
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FIG. 15: Cavity emission spectrum. The solid-thick-gray lines correspond to the situation, in which all four transitions have 
equal coupling strength, namely, a = /3 = 1, gci = Qib = g, while the dashed-thick-gray lines give the cavity emission spectrum 
in absence of coupling with the quantum dot. The solid black lines are found for (a) a = ft = 0, gci = giB = \2g, S = —0.2g, 
e = 0; (b) a = /3 = 1, g G i = l-5g, g iB = O.663, 5 = -3g, e = 0; (c) a = f) = 0, gai = 2g, g 1B = g, 5 = -2.7g, e = 0; (d) 
a = (3 — 1, gci = 2g, g\B = 0.5gr, 5 = — 4.2g, e = 0. The crosses are found by setting e — g. The other parameters are x = 20g, 
7 = 0.05(7, r<j = 0, P = 1.8(7, fc = O.lg. The effective two-photon coupling is the same for all the curves: g e « = —g 2 /x- I n 
each plot it is explicitly indicated the steady state average photon number corresponding to each curve. The insets depict the 
coupling scheme for each set of parameters. 



transition depends on the number of photons. However 
as far as the condition S 2 p + S n (aJa) <C x is fulfilled, 
so that the perturbative expansion is valid, then the de- 
tuning 6 can be set in order to minimize the effect of 
the light shift. In this case equal values of the effec- 
tive coupling (? ff give similar dynamics independently of 
the exact values of the coupling constants gci, 9G2, 9ib 
and g2B- Moreover, the finite value of e affects only the 
single exciton dynamics (see also Eq. IA3|) , but it is not 
relevant for the two photon transition, provided that the 
perturbative treatment, allowing for eliminating the in- 
termediate levels, is valid. 



In Figure [15] we report the cavity emission spectrum 
for various cases. In (a) we compare the cavity emission 
spectrum obtained with the simplified model of Sec. |TT] 
with the spectrum found when only one exciton state 
is coupled to the cavity; in (b) we study the situation in 
which the coupling with the ground-exciton transitions is 
larger than that with the exciton-biexciton transitions; in 
(c) only one exciton state is coupled to the cavity and the 
couplings to the ground and biexciton states are different; 
finally (d) is found by increasing the difference between 
the coupling with the ground-exciton and the exciton- 
biexciton transitions. In all the cases we consider both 
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the situation in which e = and 6 ^ 0; moreover the 
parameters are set in order to keep g c g equal for all the 
curves. The detuning 5 for each curve is found by maxi- 
mizing the average steady-state cavity photon number as 
a function of 5 for the corresponding set of parameters. 
We observe that similar results for the cavity field are 
found in Fig. [T5] (a), (b) and (c). On the other hand, 
a relevant difference is found in (d) when the couplings 
with ground-cxciton and exciton-bicxciton transition are 
sufficiently separate in frequency. 



In conclusion, similar two photon laser dynamics can 
be obtained, independently of the exact values of the 
coupling constants and exciton detuning, as far as the 
parameters that scales the effective two-photon dynamics 
are similar. Analogous considerations are valid also for 
the effect of the spontaneous emission rates which are 
in general different for the four electronic transition: the 
main features of the laser dynamics are independent from 
the exact values of the individual rates. 
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